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I. INTRODUCTION 

Casimir-Lifshitz forces provide an experimental window on the nature of quantum vacuum at the macroscopic 
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We report on measurements of forces acting between two conducting surfaces in a spherical-plane 
configuration in the 35 nm-1 ^m separation range. The measurements are obtained by perform- 
ing electrostatic calibrations followed by a residuals analysis after subtracting the electrostatic- 
dependent component. We find in all runs optimal fitting of the calibrations for exponents smaller 
than the one predicted by electrostatics for an ideal sphere-plane geometry. We also find that the 
external bias potential necessary to minimize the electrostatic contribution depends on the sphere- 
plane distance. In spite of these anomalies, with a proper but model-dependent subtraction of the 
O ' electrostatic contribution we have found evidence for short-distance attractive forces of magnitude 

I comparable to the expected Casimir-Lifshitz force. We finally discuss the relevance of our findings 

, in the more general context of Casimir-Lifshitz force measurements, with particular regard to the 

critical issues of the electrical and geometrical characterization of the involved surfaces. 

PACS numbers: 12.20.Fv, 03.70.+k, 04.80.Cc, ll.lO.Wx 

^ . Casimir-Lifshitz forces 0, Bp: 

qh' level [1, 0, H, S 0, B EOj HH ■ After pioneering attempts in the original parallel plate configuration , and in 
a variant of this configuration based upon a sphere and a plane [l^, the claimed accuracy of recent experiments 
ranges from 15% in the parallel plane case [13] to 0.2-5% in the sphere-plane case [H, [H, [H, [I!]. In the case 
^ , of the sphere-plane configuration, the Casimir force is typically evaluated by using the so-called proximity force 
QQ ' approximation (PFA) [20l.l2l|. introducing a theoretical uncertainty estimated to be in the 0.1 % range, and recently 

\ investigated experimentally |23 |. 
CO ■ A sphere-plane experiment has been performed, with a radius of curvature for the sphere of order of cm, at 
' relatively large distances above 1/im jT5| . In this case the largest deviation from the bare Casimir formula for an 
, idealized configuration (perfect conductors, zero temperature) is expected to be the thermal contribution to the 
7-H ■ radiation pressure on the surfaces, still below the sensitivity of the apparatus at the largest explored distances. Three 
OO , sphere-plane experiments have also been performed in a quite distinct regime, with smaller radius of curvature, of 
■ order 100 /im, and in the range below one micrometer [3, [l3j Ull- At distances less than 1 /xm the correction 
\ to the Casimir force due to finite conductivity and surface roughness of the substrates cannot be neglected at the 
• '~j ■ level of accuracy claimed in these three experiments. At least other five groups have recently performed Casimir force 
/\ ' experiments in the sphere-plane configuration [23, [H, [2^, [13, [11] . 

^ . Mastering the Casimir force at the highest level of accuracy is essential to provide reliable limits to other macro- 
■ ■ ■ ' scopic forces acting in the micrometer range, such as corrections to the Newtonian gravitational force independently 
predicted in many attempts to unify gravitation with the remaining fundamental interactions (29l . [sol . [sij . Although 
in the nanometer range the Casimir force loses its universal nature morphing into more specific, structure-dependent, 
molecular van der Waals forces, its study is sometimes considered of some relevance for designing nanomechanical 
structures and to investigate nonlinear effects, as first discussed in [s^ and then experimentally demonstrated in 

MM- 

In this paper we report more extensively on measurements performed in a sphere-plane geometry exploring a novel 
range of parameters that interpolates between the two previous sets of sphere-plane measurements [35|. We use a 
gold-coated spherical lens with large radius of curvature, similar to the experiment performed by Steve Lamoreaux 
in Seattle, while at the same time exploring distances down to few tens of nanometers from the point of contact 
between the sphere and the plane, similar to more microscopic setups using microresonators. While we do achieve 
some evidence for the expected Casimir-Lifshitz force, our measurements suggest that a reanalysis of systematic effects 
in previous experiments will be beneficial to assess the accuracy with which this force has been measured so far, with 
particular regard to the role of the residual electrostatic forces [s^ [s^l . 

The paper is organized as follows: in Section II we describe the experimental apparatus and the basics of our 
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FIG. 1: Schematic of the experimental setup and the data acquisition scheme. The silicon resonator is opposed by a spherical 
mirror, both Au coated. The latter can be moved along the vertical direction by driving two mechanical actuators for coarse 
approach, and with a piezoelectric actuator for fine control of the separation. The outside face of the resonator reflects light 
coming from a fiber optic interferometer, located 20-100 /^m above it. The resonator is driven weakly at its resonant frequency 
with a piezoelectric actuator connected in a phase-locked loop (PLL) to the optical fiber. 



measurement technique. In Section III we describe a procedure for electrostatic calibrations to determine various 
critical parameters such as the calibration factor, the absolute separation, and the minimizing potential related to the 
contact potential. Careful control over the system parameters allows us to search for a non-Coulombian contribution to 
the observed force signal and therefore to test this residual against various hypothesis such as a possible uncompensated 
electrical voltage, and the presence of the Casimir-Lifsthiz force, as described in Section IV. In Section V we critically 
assess systematic effects in our measurements and compare our findings to both previous short-range Casimir force 
measurements and long-range atomic force microscopy measurements. 



II. EXPERIMENTAL SET-UP 

The experimental setup is an upgrade of the one already described in a previous paper [sl] for measuring the 
Casimir force in the cylinder-plane configuration as first theoretically discussed in [s^. Forces acting between a flat 
surface and a spherical surface are detected by measuring the shift in the mechanical oscillation frequency of a silicon 
cantilever that acts as the plane surface. The main components in this setup are the use of a phase-locked loop scheme 
j40| to drive the cantilever at its resonant frequency, the use of optical quality surfaces for both the resonator and the 
sphere, and a faster and more flexible data acquisition scheme. 

A schematic of the apparatus is shown in Fig. 1. The rectangular cantilever, of length L = (22.56 ± 0.01) mm, 
width w = (9.93 ± 0.01) mm and thickness t = (330 ± 10) /im was laser-cut from a silicon wafer. Given the density 
of silicon, p = 2.3 x 10'^ kg/m^, the physical mass of the cantilever is rrip = (1.72 ± 0.05) X lO"'^ kg. The cantilever 
is sandwiched between two aluminum clamping structures, electrically isolated and thermally stabilized by a Peltier 
cooler. This last feature is required since we have observed that the resonant frequency of the cantilever drifts as 
much as 2 Hz as the ambient temperature fluctuates about 1 K over the duration of a typical run of a few hours. 
The effect of the thermal drifts has also been mitigated off-line by a proper implementation of the data acquisition 
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FIG. 2: Calibration of the PZT actuator governing the sphere-plane separation. Four calibrations of the PZT actuator have 
been performed during the overall data taking, of which only two are reported here for simplicity, indicated with red-square 
and black-dot points, and their best fit with a sinusoidal function (red, dotted line, and black, continuous line, respectively). 
The average value of the actuation coefficient from the four data sets is P — (87 ± 2) nm/V. 

sequence, as described in Section III. 

Below the cantilever a spherical mirror with radius of curvature R ~ (30.9±0.15) mm and diameter a — (8.00±0.25) 
mm is mounted on an aluminum frame connected to two motorized actuators allowing for coarse translational motion, 
plus an additional piezoelectric transducer for fine translational motion driven with a bias VpzT- The spherical lens is 
gold-coated by the manifacturing company with a thickness of 250 nm (± 10 %) The surface of the cantilever facing 
the spherical lens is coated with a 200 nm layer of gold (± 10 %) by thermal evaporation, with the rate of evaporation 
kept below 50 A/sec at a pressure of 3 x 10"^ mbar to ensure better conditions for homogeneous coating. 

The predicted frequency of the fundamental flexural mode of the cantilever is z/p = {0.162t/L^)y^E/ p with E the 
Young's modulus of the material [4l] |. In our case this yields a resonant frequency around 894 Hz, within few percent 
from the measured frequency at the typical vacuum pressure of 1.6 x 10~^ Torr in the vacuum chamber. The stiffness 
k of the resonator can be estimated by assuming an effective mass for the resonator mode roughly equal to the physical 
mass, which yields k = 1.036Ewt^ / ~ 5.4 x 10'^ N/m. The stiffness is three or four orders of magnitude higher than 
the typical cantilevers used in atomic force microscopy, implying both advantages and drawbacks. A large stiffness 
allows us to reliably achieve small gaps with little static deflection of the cantilever in the presence of strong electric 
fields for calibration purposes, at the expense of a lower force sensitivity. By decreasing the stiffness of the cantilever 
both the reachable minimum gap before snapping between the two surfaces occurs and the largest explorable gap for 
a given signal to noise ratio are increased. 

The motion of the cantilever is detected by using a fiber optic interferometer [4^ 1 positioned a few microns above 
the upper cantilever face and fed by temperature-stabilized diode laser with an adjustable power in the 5-10 mW 
range at the wavelength of 781 nm. The output signal from the interferometer is filtered and amplified through a 
single reference mode lock-in amplifier and is fed back into the piezoelectric actuator driving the cantilever motion. 
The phase between the input and the output signals is properly chosen to maximize the vibration amplitude, typically 
around 30-40 degrees. This scheme is much faster and more efficient than the previous open-loop scheme as described 
in [111 in which white noise was used as the driving source without a feedback, and the complete FFT spectrum was 
acquired. The vibration frequency was measured by a counter with 10 mHz resolution, 100 times faster than in the 
open-loop scheme. The front side of the chamber has a large viewport allowing for visual inspection of the relevant 
components through an optical microscope, as shown in Fig. 10 of This also provides for a quick assessment 
of the optical fiber location, which can be manually adjusted within few hundred micrometers from the underlying 
resonator using a feedthrough micrometer. 

The PZT was calibrated by using the same fiber optic interferometer with the cantilever removed and the spherical 
mirror replaced by a flat mirror. Assuming a linear relationship Ax = /3AVpzt, the DC interference amplitude versus 
the voltage applied to the PZT, as shown in Fig. 2, was fit with a sinusoidal function and the distance inferred from 
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its period to be half of the 781 nm laser wavelength. The average of four separate calibration runs yields (3 = (87 ± 2) 
nni/V, and it was found to be independent of the voltage applied to the PZT, within the error, in the 0-100 V range. 
The sensitivity of the fiber optic interferometer is evaluated in the same configuration by driving the PZT with a 
monoschomatic signal at various amplitudes, and looking the corresponding signal at the FFT spectrum analyzer. 
When the driving amplitude at a frequency of 1 KHz (still well below the maximum response frequency of the PZT 
of 69 KHz) is reduced at the level of a peak-to-peak signal of 0.25 mV, the measured amplitude at the FFT spectral 
analyzer becomes equal to the FFT broadband noise level of 10 /iV^ms with an integration time of 10 s. Using the 
determined value of the actuation coefficient of the PZT, we estimate the minimum detectable signal with SNR ~ 1 
to be ~ 0.80A/\/Hz, similar to the one reported in [i^ . 

For a generic distance-dependent force and for a small amplitude of the cantilever oscillation (2-3 nm as inferred 
from the interferometer signal) with respect to the average separation between the two surfaces x, the measured 
frequency of the resonator z^,„ is shifted with respect to the proper frequency of the resonator Vp in such a way that 

where nics is the effective mass of the mode of oscillation of the resonator. In general, the total force acting on the 
resonator is the aggregate from independent sources and the square of the measured frequency Vm of a cantilever can 
be identified as the sum of contributions from those sources, 

i.^(x, V) = + ^iJlix, V) + Avlix), (2) 

where A.v'^ is the frequency shift due to electrical force gradients for instance due to an external bias voltage V , 
and A^-^ is the frequency shift subject to force gradients of non-electrostatic nature, for instance the Casimir force, 
or random drifts of instrumental and environmental origin. In the next section we will discuss the response of the 
resonator in the case of electrostatic forces. 

III. ELECTROSTATIC CALIBRATIONS 

We start with the electrostatic energy stored in a generic capacitance C{x) biased by an external potential V , and 
under the presence of a contact potential Vc{x) which in general may depend on the distance: 

E,i^^[V-V,{x)f. (3) 

As recently pointed out by Lamoreaux [i^ , if the contact potential depends on distance the electrostatic force will 
include an additional term independent upon the external bias voltage. Since in our configuration the observable is 
the frequency shift, proportional to the second derivative of the electrostatic energy with respect to x, we have, after 
regrouping the various terms in powers of 1/ — 14: 

Ki -^^{V~V,f^ [2CX. + {V^V,) + CVi^ . (4) 

This can be regrouped again as i?"; = A{y — Vc + B)^ + D and, identifying the coefficients A, B, D to match the 
terms in Eq. 4, we obtain 

This allows to write: 

^o(x) = F.(x) + ^^^^i±^ (7) 

is the minimizing external potential, i. e. the external bias minimizing the electrostatic contribution to the frequency 
shift. It is clear that, even if the external bias voltage is chosen to minimize the electrostatic contribution, an 
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FIG. 3: Parabolic dependence of the squared frequency vl^ upon the external bias voltage V for a typical experimental run. 
The different curves correspond to different values of the sphere-plane separation (defined in the legend in terms of the voltage 
applied to the PZT used for approaching the sphere). The range of external bias voltages is adapted for each sphere-plane 
separation to maintain comparable values of the maximum detected frequency shift, of about 1 Hz. This corresponds to keep a 
constant value for the electric field in the gap. The best fit with a parabolic function is also shown, with the parabolic curves 
relative to the largest PZT voltages (corresponding to smaller gap separations) displaying largest curvatures. No significant 
dependence of the fitting parameters has been observed by changing the external bias voltage span. 



irreducible term containing the first and second spatial derivatives of Vc in Eq. 7 still causes a frequency shift. The 
expected frequency shift due to Coulombian interactions is therefore: 
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(8) 



We identify three contributions on the RHS of Eq. 8: the first term is unrelated to external or internal voltages, 
the second term depends on the external voltage, and the last term is a Coulombian contribution only related to the 
spatial dependence of the contact potential. The electrostatic calibrations consist in the study of the response of the 
cantilever frequency to an external potential, therefore selecting the second contribution alone which allows to infer 
both the effective mass of the resonator and the minimizing potential Vq. This knowledge is then used to infer the 
contact potential assuming a specific form of the functional dependence on distance of Vq, which is in turn used to 
identify the third contribution in Eq. 8. The subtraction of the total Coulombian contribution allows then to evaluate 
the square frequency related to non-Coulombian forces and to random drifts due to background noise in the apparatus 
(see section IV) as first discussed in AFM measurements by and in the specific context of Casimir measurements 
in [H. 

Regarding the capacitance, for the sphere-plane configuration and for our choice of parameters, the proximity force 
approximation (PFA) holds and the capacitance is C{x) = 27reoi?ln(i?/x). This formula in principle must be corrected 
for the finite spherical mirror diameter used in alternative to a full sphere, but this is a sub-leading PFA correction 
approximately equal to 0.1% for a typical separation of d — Ifim. Therefore, the measured frequency Vm of the 
cantilever can be parameterized as 



iylix)^iyU^.)-IU{x)iV-yof 



(9) 



where Vq{x) = vi + ISxr^ix, Vq) + Ai/^(x), a parabola whose maximum is reached when the applied voltage at a given 
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FIG. 4: Electrostatic curvature coefficient versus the voltage applied to the PZT actuator for four different runs, and related fits 
with the expected fixed exponent e = —2 (blue dashed line) as in Eq. 10 and by allowing the exponent to be a free parameter 
(red continuous line) . In run 1 the location of the surface of the sphere in proximity of the cantilever is different from the other 
runs, as the latter were taken after tilting the sphere by an angle equal to 5xl0~^ radians. 



distance equals to the minimizing potential Vo. The parabola curvature K(.i{x) = eaR/ATrrricgx'^ reflects the cantilever 
response to externally applied electric forces at a given distance x and allows to extract the effective mass of the 
resonator mode, and it can be measured by sweeping the applied bias over a bipolar range to discern the parabolic 
dependence of the frequency shift. The minimizing potential Vq is found by looking at the horizontal displacement 
of the parabola. In Fig. 3 we show examples of electrostatic calibrations in which the square resonator frequency 
is plotted versus the bias voltage for different values of the gap separation between the sphere and the plane. In 
order to minimize the resonator frequency drifts due to fluctuations of ambient temperature, we measure the unbiased 
frequency before and after each measurement taken at a given bias voltage V. The average of the two unbiased 
frequencies is then used as the reference to evaluate the shift. We notice two remarkable features in Fig. 3. First, 
the peak value of the frequency is different in the four cases, reflecting the possible presence of a distance-dependent 
force combined with possible drifts in the value of the intrinsic frequency. Second, the location of the peak value 
of the resonator frequency in the parabolic fit does not occur at the same value of external voltage, indicating a 
distance-dependent minimizing potential. 

The absolute distance between the sphere and the plane is not known with sufficient accuracy prior to the force 
measurement. In the experiment, the gap is varied by the voltage applied to the PZT (Vpzx) and consequently Ke\ 
is a function of the relative distance (i.e., of the applied Vpzx)- This requires an additional fitting parameter Vp2Xj 
which would cause contact and must be inferred from fitting the function 
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where a is a calibration factor containing the effective mass of the cantilever through a = eoR/4:TTmcs(3'^ . In Fig. 4 
the curvature parameter K^i is plotted versus the piezoelectric actuator voltage T^pzT for four different experimental 
runs. Assuming the PZT actuation is linear, the absolute separation for a given VpzT can be inferred in two ways: 
a;(VpzT) = PiV^zT - Vpzt) or a;(VpzT) = /3[Q;/i^ci(VpzT)]^/^- Therefore, the absolute distance can be inferred either 
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FIG. 5; Plot of x{Vpzt) converted into distance, as evaluated from the calibration factor a for the inverse square law, versus 
the same function evaluated from the asymptotic limit of the fit function. The left plot uses the expected electrostatic law with 
an exponent e — —2, while the right plot is obtained by optimizing the exponent in the curvature scaling with distance. Red 
lines are added to provide an eyeguide for the deviations from the expected consistency between the two determinations of the 
same quantity. 

from the asymptotic limit Vp'2x of the fit function or from the calibration factor a of the same function, indicating an 
intcrdependency of the two physical parameters appearing in Eq. (|10[) . Note that quite large values of Kd, of order 
—2 X in runs 1. 3, and 4, have been achieved. These large signals serve to accurately constrain the fit 

parameters owing to the high stiffness of the cantilever, which allows us to explore very small distances between the 
sphere and the cantilever without stiction. 

Comparison of the two methods for determining x^Vpzt) in Fig. 5 (left plot) reveals disagreement when applying 
the inverse square law in Eq. 1101 Direct measurement of the sphere-plane capacitance versus gap provides a third 
assessment of Vp2T- In Fig. 6 the capacitance versus PZT voltage is shown with the fitting curve corresponding to 
the capacitance between a sphere and a plane in the PFA approximation. 

Another issue complicating the electrostatic calibrations is evident in Fig. 7 displaying Vq inferred from Fig. 3 
versus Vpzx- While the error bars are relatively large at large gap separations, they become smaller at smaller gaps, 
and two distinct dependencies are observed. In runs 1 and 2 the minimizing potential seems to depend linearly on 
the gap separation at small separations, while in the remaining runs it appears to tend towards a constant value. 
The significance of this dependence in the context of high-precision measurements of the Casimir force is discussed in 
section V. 

The error in the curvature parameter K^i has been determined with dedicated runs in which conditions nominally 
leading to the same gap size are maintained, and the result is shown in Fig. 8. The overall time drift in K^i, if converted 
into a thermally induced change in gap size, corresponds to ±200 nm. The average random uncertainty in Kgi has 
been evaluated by taking the data of Fig. 8 and subtracting for each point the long term drift by means of a moving 
average with a defined time window. By choosing the average time window equal to 4 curvature measurements (each 
measurement taking 8-10 minutes depending on the number of points obtained in the frequency versus bias voltage 
measurement) we obtain a relative error in K^i equal to about 4% error on top of the fitting uncertainty associated 
with the parabolic fit. This allows us to quantitatively compare fitting functions for the distance dependence of Kd, 
and has resulted in another anomalous behavior. Indeed, the reduced is near one when the exponent is fit in the 
-1.7 to -1.8 range, while it is at least one order of magnitude higher for fixed -2.0 exponent, as shown in Table 1. 
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FIG. 6: Capacitance versus PZT voltage data (black circles) and its best fit (red line). The fit is based on the PFA formula 
C{d) = Co + A ln[f3{VpzT ~ Vpzt)], where the actuation coefficient (3 = (87 ±2) nm/V. The best fit yields, for the coefficients Co 
and A, the values Co ~ (193.9±0.2) pF and A — — (1.757±0.002) pF, this last in agreement within two standard deviations with 
the less accurate theoretical expectation A — —2TieoR = —(1.72 ± 0.02) pF. The same fit provides a value for the asymptotic 
value of the PZT voltage corresponding to zero gap distance the value Vpzx = (69.31 ± 0.01) Vm, equivalent to a distance of 
47 ± 2 nm. The best fit corresponds to a reduced — 2.9. 



Adherence of our data to a strict power law from the farthest to the closest approach in all four cases implies the drift 
played little role during our data runs, except possibly in the case of Run 4. However, we find that the data clearly 
fail to follow the inverse square law dependence of the electrostatic coefficient upon the sphere-plane separation. More 
quantitatively, if the fitting exponent is left as a free parameter, our experimental data from four separate sequences 
follow a power law with exponents — 1.70 ± 0.01, —1.77 ± 0.02, —1.80 ± 0.01, — 1.54 ± 0.02, far from the expected value 
of -2. When substituted into Eq. [TOl these exponent values produce better agreement between the two methods for 
determining a;(VpzT), as shown in the right plot in Fig. 5. 

The stability of the electrostatic result has been checked through repeating the data fit of Ke\ versus VpzT, starting 
with few points at the largest distances and by progressively including the data point corresponding to the closer 





Run 1 


Run 2 


Run 3 


Run 4 


fixed exponent 


-2 


-2 


-2 


-2 


Q (Hz^V-^m^) 


6200 ± 98 


6197 ±97 


6701 ± 101 


6438 ± 156 


xo (nm) 


64.4 ± 1.7 


90.5 ±2.6 


62.6 ± 1.7 


93.1 ±2.6 


xVdof 


15.9 


7.7 


6.9 


36.5 


free exponent 


-1.70 ±0.01 


-1.77 ±0.02 


-1.80 ±0.01 


-1.54 ±0.02 


a (Hz^y-^m'') 


2805 ± 92 


3021 ± 153 


3732 ± 144 


1415 ± 69 


Xo (nm) 


29.6 ±0.9 


49.6 ±1.7 


35.7 ± 1.7 


20.0 ± 1.7 


xVdof 


1.0 


0.8 


1.2 


7.0 



TABLE I: Fitting parameters of the electrostatic calibrations, and values of the reduced , as the absolute normalized 
by the number of degrees of freedom (DOF) in the fitting for a fixed exponent equal to 2 as expected from electrostatic 
considerations, and for a free exponent e. Note that the calibration factor a characterizes the sensitivity of our apparatus, 
yielding the effective mass in the case of fitting with the expected Coulomb exponent equal to 2. Variations of a from run to 
run are attributable to the rearrangement of the resonator clamping system, which strongly affects its effective mass. If the 
exponent is different from the expected value of 2, then the dimension of a changes and the effective mass cannot be properly 
evaluated. 



9 







-0.1 

> -0.2 

> -0.3 
-0.4 
-0.5 



"I I r 



A.I 




Run 3 



J I I I I I I I I I L_ 



10 20 30 40 50 60 
Vp2-p(V) 




10 20 30 40 50 60 



Vp2.-p(V) 



FIG. 7: Plot of the minimizing potential versus Vpzx for the same four runs shown in Fig. 4. Two behaviors are evident at 
small separation (large Vpzt)'- either an approximately constant potential (Run 3 and Run 4) or a potential linearly dependent 
upon the separation (Run 1 and, to a less extent, Run 2). Note that in run 1 Vb is positive while for Runs 2-4 it is negative 
and converges to Vb — -150 mV at the smallest separations. 



distances. In this way the intrinsic instabihty of the fit will manifest itself through a systematic variation in the fit 
parameters subject to a choice of data points included. Figures 9 and 10 show the result of the stability test for Run 
1-4, revealing dependence of the fits on the number of points considered. As we will discuss later in Section IV, this is 
limiting the accuracy of the determination of the Casimir force. In Table 1 we report the values of the parameters of 
the best fit of the electrostatic calibration for the four runs, leaving the exponent free, and with the exponent frozen 
at the expected value of -2. In the latter case one can infer the effective mass of the cantilever and this turns out to 
be larger by an order of magnitude than the physical mass, against the expectation to be a fraction of the physical 
mass [4^, [43, EBl- Apart from the uncertainties related to the clamping of the resonator, this is another signal that 
the calibration fitting with a pure electrostatic contribution neglects some systematic effect. 

The unexpected power law poses a significant limit on the validity of our electrostatic calibration. Among the 
possible systematic effects causing this deviation from the expected Coulombian behaviour, we have considered the 
following. 



A. Static deflection of cantilever 



The spring constant of our cantilever is extremely stiff (about 5400 N/m). Using Hooke's law, a deflection experi- 
enced by the cantilever due to an electrostatic force at 100 nm with an applied voltage of 100 mV is less than 0.2 A. 
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FIG. 8: Long-term drift of the curvature coefficient Jfci during a 12 hours test. The vertical bar represents the equivalent drift 
in terms of the change in the sphere-plane separation, in the conservative hypothesis that all the drift is attributable to the gap 
drift. Subtracting each point with a window average of 4 measurements one gets a relative error in K^i of 4 %. The sensitivity 
to the choice of the average window time has been also assessed, with relative errors of 2.9 %, 4.1 %, 4.4 % and 4.6 % for 
alternative choices of 2, 6, 8, and 10 measurements, respectively. 



Hence, the static deflection should play Uttlc role. 



B. Thermal drift 



Even though the temperature of the cantilever is actively stabilized by a Peltier cooler, the rest of the system is still 
subject to global thermal variation. In order to see this, we have measured fcei with respect to time at a nominally 
fixed distance. In the worst circumstance, the gap separation during the course of measurements can drift as much 
as 200 nm in either direction. Although such a drift could in principle affect the inferred exponent, a highly unlikely 
non-linear monotonic drift would be necessary to account for the consistently observed anomaly in each independent 
run. 



C. Nonlinearity of the PZT actuation 



The linearity of the PZT translation has been tested under a number of different circumstances. Notice that the 
translation intervals between the data points in each of the runs shown in are completely random. Yet, all of the runs 
obey a specific power law in all distances. The PZT was also independently calibrated by means of the fiber optic 
interferometer with a consistent, linear actuation coefficient factor /3 = 87± 2 nm/V. Also, the presence of significant 
nonlinearities should manifest also as deviations from the PFA fitting of the sphere-plane capacitance versus the PZT 
voltage. As visible in Fig. 6, this class of measurements strongly constraint the PZT nonlinearity. 
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FIG. 9: Test of stability of the fit of the electrostatic calibrations: shown is the determination from the fitting of the calibration 
factor a in the case of the fixed -2 exponent versus the sphere-plane distance for the fixed -2 exponent (in units of Hz^ 
m^) expected from electrostatics (diamond points, black), and for the best exponent left as a free parameter (circle points, 
red). The number of points used in the data fit includes up to the gap distance evaluated at that point. The best exponent 
points show convergence to a constant value at the closest distances in all four runs, unlike the electrostatic case with the fixed 
exponent. 



D. Nonlinear oscillation of cantilever 

The cantilever is driven at resonance in a phase-locked loop, a routine technique adopted by many groups [isl . 
HEj l49l . [sot . Higher order terms in the force expansion should produce higher harmonics of the drive frequency. 
Then the assumption that the frequency shift is simply proportional to the gradient of the external force could lead 
to erroneous assessments of /cci, eventually affecting the exponent. We have not observed higher harmonics in the 
frequency spectrum of the resonator. 



E. Deviation from geometrical ideality 

The deviation from geometrical ideality and its influence on the local capacitances [sij could in principle play a 
role, both in having a different shape or a variable radius of curvature, and in the deviations from the PFA formula. 
In our case, even at the smallest explored distances, the PFA correction due to the roughness is estimated to be of 
0.4% in the electrostatic calibration, well below our sensitivity. A more careful analysis beyond the PFA is possible 
taking into account the actual surface profiles. It is however worth to point out that at least the simple hypothesis of 
a non uniform radius of curvature of the spherical surface cannot take into account the anomalous scaling. In PFA, 
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FIG. 10: Test of stability of the fit of the electrostatic calibrations: shown is the determination from the fitting of the offset 
distance xq ~ PVp^T versus the sphere- plane distance in the case of the fixed -2 exponent expected from electrostatics (diamond 
points, black), and for the best exponent left as a free parameter (circle points, red). Clearly, the trend associated with the 
fixed exponent displays a wide variation in the fitting parameter Vp^Tj while the overall behavior driven by the best exponent 
shows a degree of stability in all distances. 



the force signal (for both electrostatic and Casimir) is directly proportional to the radius of curvature and the total 
force will be therefore obtained as an integral along the radii of curvature with no deformations to the inverse square 
law for the distance dependence. In [s^l an interesting geometry has been discussed which could take into account 
the anomalous exponent we have observed, however the stringent agreement between the capacitance measurements 
and the PFA expression for a pure sphere-plane geometry rules out this hypothesis. 



F. Surface roughness 



Roughness corrections become important at small separations between the surfaces, both for electrostatic cali- 
brations and for analyzing their residuals in the search for Casimir forces. Usually roughness corrections in the 
electrostatic calibration are disregarded as such calibrations are performed at large distances. Here we consider the 
simplest analysis based on the proximity force approximation. Let us assume that both surfaces have stochastic 
roughness with rms amplitudes (/i^) and (hp) for the sphere and plane respectively. Apart from the usual PFA 
condition to treat the curvature effects due to the spherical lens (c? <JC R), we further assume that the sphere-plane 
distance d is much smaller than the lateral roughness correlation length ^ on each surface (this is the condition for 
the applicability of PFA to roughness considerations, only valid for very smooth surfaces). Further assuming that the 
rms roughness amplitudes are the smallest lengthscales in the problem, the PFA second-order perturbation correction 
to the sphere-plane electrostatic force is 
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FIG. 11: Atomic force microscopy images of the cantilever surfaces, of size 10 x 10 /im (left) and 1x1 fim (right). The overall 
rms roughness for the larger surface is 1.54 nm (with a peak in the profile due to dust of 4.78 nm) and for the smaller surface 
is 0.79 nm (with a peak in the profile of 2.2 nm). 



where Fq is the same force evaluated for smooth surfaces. Figure 11 shovirs typical images of the cantilever surface 
taken using AFM microscopy. We measured the surface profile for our surfaces and estimate the rms roughness of 
{h1)= 4 nm^ and {hp)= 2.4 nm^. The resulting deviation from ideality is below the estimated sensitivity of the 
apparatus, and an attempt to fit the electrostatic calibrations including the roughness at the PFA level as in Eq. 
[TT] did result in significantly larger reduced ■ Roughness corrections to the electrostatic calibration potentially 
represent a systematic effect that could influence the residuals in the search of the Casimir contribution. 



Patch effects are expected to induce deviations from the Coulombian scaling with distance [37| , and the fact that we 
have found in some calibrations a distance dependent minimizing potential could validate this hypothesis. However, 
this docs not explain yet while the anomalous exponent is systematically observed in all runs, while only in some of 
them a manifest distance-dependent minimizing potential occurs. 

IV. ANALYSIS OF THE RESIDUALS OF THE ELECTROSTATIC CALIBRATIONS AND TESTS OF 

THE CASIMIR FORCE 

In the previous section we have found that the electrostatic calibrations performed in the entire measurement 
range show the presence of a distance-dependent minimizing potential (and therefore a distance-dependent contact 
potential) and an apparent deviation from the expectation for the sphere-plane capacitance dependence on distance. 
While the first systematic effect can be taken into account in the data analysis, the anomalous distance scaling of 
the electric force and the lack of its complete characterization prevents in principle to master the calibrations. This 
prevents the reliable determination of the relevant parameters of the system, for instance its effective mass, and 




(11) 



G. Patch potentials 
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therefore the derivation of the residuals containing the Casimir force signal. In the absence of a complete analysis, 
and waiting for further experimental and theoretical input on the anomalous exponent observed in our configuration, 
in this section we analyze the cantilever square frequency f g (x) obtained by subtracting to the electrostatic calibration 
data the Coulombian electrostatic contribution based on the fitting procedure, as first performed in [4^ . This off- 
line procedure is necessary as we have evidenced a dependence on distance of Vq in various runs, and therefore we 
cannot assume a priori a unique counterbias voltage for an on-line subtraction of the residual electrostatic force as 
typically performed in various experiments. A distance-dependent residual frequency shift implies the presence of a 
non-Coulombian force or drifts in the intrinsic resonator frequency or a combination of both. It is evident in Eq. 
8 that the electrostatic residual, Az/g (x, Vb), is zero only if the contact potential is constant. This term must be 
evaluated by first inferring the contact potential Vc through numerical integration of the minimum condition Vq. 
Based on Eq. [7] we have a second order differential equation for Vc{x): 



C d'^VJx) dVc ,,,,,,, , 

where Vq{x) is determined experimentally. For the sphere-plane case in the PFA this results in the following differential 
equation: 



l^^R/^)f^ _ 2x^^ + V,{x) = V,{x). (13) 
ax ax 

to be solved numerically for instance using a Runge-Kutta integration method, once the experimental data on Vq{x) 
are interpolated with an analytic function. In the absence of a model for the dependence on distance of Vb , different 
equivalent empirical functions can be used to fit the actual data. This obviously is a strong limitation to the predictive 
power on the residuals, but at a phenomenological level has to be considered as a conservative and realistic procedure 
to infer some physics beyond the Coulombian contribution. We use, as first guess function, an exponential one such 
as Vo{x) = Vq + Ay[l — exp(— x/A)] with Vq representing the potential when the two surfaces are in contact, Vq + Al/ 
is the asymptotic value of the minimizing potential at large distances, and A the characteristic lengthscale on which 
the minimizing potential varies. Alternatively, we may use the function Vq{x) = T4og -|- AX^og ln(a;/A), characterized 
by similar parameters Mog, AViog and A, although in this case a defined asymptotic value is not available. In Fig. [T^] 
we show the sensitivity of Vc obtained through integration of Eq. [13] to the choice of the function interpolating the 
Vq distance dependence. 

Another issue requiring careful attention regards the boundary conditions chosen to solve the second order differ- 
ential equation for Vc{x). At large distances Vc{x) is expected to be constant and thus Vc{x) — Vq{x). Also, the first 
derivative of Vc{x) should be zero in order to satisfy the differential equation. Therefore we choose the boundary 
condition Vc(x„) = Vo(a;„); V^{xn) = 0, where Xn is the largest distance measured in the data set. In the case of the 
exponential fitting of Vb(x), the curve becomes quickly fiat as x increases and therefore the use of Xn at very large 
distances is a good approximation. For the logarithmic fitting this is more problematic, making the solution for Vc{x) 
quite sensitive to the choice of a;„. 

Based on this input, one can evaluate now the electrostatic contribution independent on the external bias potential, 
Aj/g(x,Vo), and further subtract it from electrostatic contribution dependent on the external bias evaluate at the 
minimizing potential, I'mi.^^ ^ — ^) {i-^-, the maximum of each parabola in Fig. 3). The result is shown in Fig. 13, 
where the dashed and dotted-dashed curves are relative to the choice of the Vq dependence on distance, as given in 
Fig. 12 (either exponential or logarithmic). We see that the plotted quantities v'^^{x,V = Vb) — /S.v1{x,Vo) indicate 
that there are further distance-dependent residuals + Azy^(a;) of non-electrostatic nature. In this analysis it is 
critical to have the most accurate control on the dependence of Vb on distance. In case of runs 2, 3, and 4 such a 
control is quite limited, and neither an exponential function nor a logarithmic function adequately describe its spatial 
dependence. In the remainder of this section we will concentrate on run 1. 

Given the serious limitations of our electrostatic calibration, it is difficult to perform a rigorousanalysis of these 
extra residuals, as they could be Casimir-Lifshitz forces, patch potential forces, etc. In the following we will assume 
that these extra residuals are only due to the Casimir-Lifshitz force, and evaluate the corresponding prediction for 
i/p -|- Ai/Qjjg(a;). As a first attempt, one can fit either the exponential (dashed line) or the logarithmic (dotted-dashed 
line) residuals in Fig. [T3]with the plain Casimir formula for zero temperature and perfect metals. In the proximity 
force approximation, the sphere-plane Casimir force is -Fcas(a;) = 27ri?£'pp (x) , where £'pp(a;) is the Casimir energy 
per unit area for the plane-plane configuration, E^p(x) = —hcir^ /720x^ . Therefore, the frequency shift due to the 
Casimir force in the sphere-plane case is Ai''^^^{x) ~ — A'cas/a;"*, with A'cas ~ nhcR/ASOmeS- Fitting the above curves 
with 1/2 - i4:cas/x^ we obtain = (790440 ± 19.3)Hz^ ATcas = (5.5164 ± 0.3390) x IQ-^^HzVm^ for the dashed 
line (exponential formula for Vb(x)), and = (790450 ± 18.3)Hz2, i^cas = (6.5795 ± 0.3220) x lO^^^Hz^m^ for 
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FIG. 12: Data for the minimizing potential Vo versus sphere-plane distance for run 1. The dashed curves represent the 
best fits of Vo assuming an exponential dependence (red) and a logarithmic dependence (black). The best fit parameters 
are Vo = (0.011 ± 0.007)V, AV = (0.25 ± 0.01)V, and A = (703 ± 93) nm in the case of the exponential function, and 
Viog = (0.07 ± 1140) V, AViog = (0.058 ± 0.003) V, and A = (140.4 ± 2.8) nm in the case of the logarithmic function. The 
continuous curves represent the corresponding plots of the contact potential Vc after solving the differential equation Eq. [13] 
as described in the text. 



the dashed-dotted line (logarithmic formula for Vo(a;)). Using the values for the radius of curvature of the sphere of 
R = 30.9mm, and an effective mass of nicff = 0.46g, the obtained coefficients i^cas are lower than theideal Casimir 
prediction (-ftTcas.th = 1.3 x 10~^^Hz^/m^)by about a factor 2, in line with the expectation that conductivity corrections 
could play a role. 

In order to analyze possible conductivity (and temperature) corrections to the Casimir force we now use the Lifshitz 
formalism [ssj . written in terms of the frequency-dependent reflection coefficients of the two gold surfaces. The plane- 
plane free energy is 

EMx) = ^Y.T. dyy' log(l - rl e'^y). (14) 

p m— "^"^ 

Here p denotes the two possible polarizations (TE and TM), 7 = 2TTkBTx/hc, and the prime on the summation sign 
indicates that the m = term is counted with half weight. The reflection amplitudes are given by the usual Frcsncl 
formulas, 

T-TM(?/,4m) = j_ Cf \ ' fTE[y,Cm)^ , (15) 

1 (^y^C^m )Pm \ Pin 



where = 2T:kBTm/h are the Matsubara frequencies, = y/rwy, and s,„ = \/ e{i£,m) — 1 + Pm- The dielectric 
permittivity at imaginary frequencies is evaluated using Kramers-Kronig relations with gold optical data from Palik 
[54| . extrapolated to low frequencies using a Drude model with parameters: iUp = 7.5eV for the plasma frequency, and 
7p = 0.061eV for the plasma relaxation parameter. The resulting theoretical prediction (using once again PFA for the 
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FIG. 13: Plot of the electrostatic residuals versus distance. The dots correspond to the maximum value of each parabola 
in Fig. 3 as a function of distance, v^{x,V = Vb). The dashed (exponential fit for Vo{x)) and dotted-dashed (logarithmic 
fit for Vo{x)) lines result from substracting from u^{x,V — Vo) the extra term Av^{x,Vo) due to the distance-dependence of 
the contact potential Vc{x). The dotted line is the Lifshitz prediction for Up + Auq^^{x). Parameters are an effective mass 
meS =0.46 g and a cantilever free oscillation frequency Vp = 889.09Hz. Optical data of gold are from Palik [sj, and Drude 
parameters are Up = 7.5eV and7 = O.OGleV. Temperature is set to T = 300K. 



sphere-plane case) for i^p + Ai''^^^{x) at T ~ 300K is plotted as the dotted line in Fig. 13, using as the free cantilever 
square frequency i/p = 7.9047767 x 10^ Hz^. It is apparent from the figure that the Casimir curve I'p + Ai'^^g(x) 
approximately follows the behavior of the curves ^'^(a;, V = Vq) — Az/g(x, Vb) (that is, the electrostatic residuals after 
substraction of the distance-dependent contact potential term). Note that Az/g(x,Vb) > 0, that means a positive 
frequency shift associated to a repulsive electrostatic residual (V^-independent) force. Had we not taken this term into 
account, the Lifshitz theory curve would noticeably depart from the electrostatic residuals v^^{x,V = Vq). 

Finally, it is important to emphasize once again that, due to the anomalous exponentwe found in our electrostatic 
calibration, we cannot make any precision measurement claims on the Casimir force after the residual analysis. Nor 
does our data have enough sensitivity to detect corrections due to temperature effects or sample dependence. On 
the other hand, it is quite possible that part of the residuals are due not only to Casimir-like forcesbut also to other 
external voltage-independent forces, not taken into account in a simple Coulomb law, such as patch potentials. 



V. SYSTEMATIC EFFECTS AND RELATIONSHIP TO PREVIOUS CASIMIR AND AFM 

EXPERIMENTS 



Our experiment spans between two different classes of experiments. In previous short-distance Casimir force 
experiments, the minimum explored distance is on the order of few hundred nanometers. This is due to the fact 
that typically, in order to increase the sensitivity at larger distances, the stiffness of the resonator is designed to be 
low. While this yields a stronger signal at large distances with respect to harder resonators, this also results in an 
earlier snapping of the resonator to the attracting surface. On the other hand, typical AFM experiments aimed at 
mapping the profile of a surface at subnanometer resolution (both depth and lateral) do not need to reach distances 
larger than 10-20 nm. We discuss here more in general the specific issues we have found that could affect at least in 
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principle the analysis of previous measurements, or at least the assessment of their accuracy. 

A. Distance dependence of the contact potential 

Our finding of the dependence of the contact potential Vc upon the distance calls for a more careful study of previous 
experiments in which the contact potential was assessed at relatively larger gaps than the one used for looking at 
the Casimir force. In order to seek for similar regularities in previous Casimir force experiments, we have requested 
electrostatic calibration data kindly provided by various groups. The result of the analysis for the data of the Lucent 



Laboratories group |17| working with spheres of much smaller radius of curvature (100 fim) is reported in the top-right 
plot in Fig. 14, showing a small but clearly visible linear correlation between Vb and the sphere-plane distance. The 
data of the AFM group in Grenoble (with microspheres of radius 20 /xm) in Fig. 14 bottom-left show also a strong 
dependence of Vq at large distances, with a more complex behavior as evident from the change in polarity of the 
minimizing potential around 5 fiui. In the case of the experiment performed at Indiana University-Purdue University 
Indianapolis (lUPUI), with a microsphere of radius (148.7 ± 0.2) fim, sec bottom-right plot, no evident dependence on 
the distance can be assessed, and in that experiment a constant counterbias equal to the average value was assumed. 
A linear trend, with a systematic difference in the contact potential in the explored range of about 3 mV around 
random fluctuations of comparable size, seems also present in Fig. 4 of (sBj using microspheres of radius 200 /im, 
although the coarse scale used for the vertical units makes harder a more careful and quantitative assessment. 

Although it would be necessary to collect more data, one noticeable emerging pattern from the comparison of 
the experiments in Fig. 14 is that experiments with intermediate radii of curvature generally seem to provide a 
milder dependence of Vb within a comparable distance range. It is plausible that spheres with very small radius of 
curvature as in the example from the Grenoble group are more vulnerable to local defects or impurities, while large 
radii of curvature (and therefore large active interaction surfaces) like in our case are sensitive to large scale variations 
in the patch potentials. It will be interesting to investigate, beside the obvious material composition, the geometry 
dependence of minimizing and contact potentials. An interesting study of the minimizing potential and its dependence 
upon distance and drift in time has been recently reported for the more macroscopic setup of a torsional balance [56| . 

To assess the impact of the distance-dependent minimizing potential on the precision of the Casimir force mea- 
surements, following [sTj . we evaluate the equivalent voltage necessary to mimic the ideal (zero temperature, perfect 
conductors) Casimir force at a given distance, and in the sphere-plane case we have 

^-(^)-120U) ^' 

implying that the Casimir force can be simulated, at a distance of 1 /zm, by just having 17.5 mV of uncompensated 
voltage between the two surfaces. For real materials the Casimir force is weaker, and can therefore be mimicked by 
an even smaller uncompensated voltage. If the minimizing potential Vq is independent of the distance over the entire 
range of explored distances, it is legitimate to use an external counterbias to cancel its effect. Otherwise, a fixed 
counterbias will only cancel out the minimizing potential at a given distance and, depending on the slope of Vq versus 
distance at smaller gaps, most if not all of the measured force could be due to the uncompensated potential. Some of 
the determinations of the contact potential in previous experiments have been performed at relatively large gaps with 
respect to those at which the Casimir forces observations are reported. In some cases, in particular for low stiffness 
microcantilcvers, this is necessary because electrostatic forces at small distances will cause snapping or instabilities, 
unless very small bias voltages, limited by the voltage supply specifics, arc applied. This suggests that a reanalysis 
of the data collected on Casimir forces so far will be beneficial in the light of these findings, also including a more 
careful scrutiny in the next round of measurements. Notice also that to have for instance a 0.1 % precision in the 
determination of the Casimir force the electrostatic background should be controlled with an accuracy better than 
550 /iV at 1 fim, which it is also equivalent to perform precision electrometry. In other words, the precision of 0.1 % 
is equivalent to be able to discern a minimum voltage of 550 fiY applied to the gap, a very difficult task considering 
the presence of comparable voltage drifts and patch potentials, although other groups have instead concluded that 
the estimated patch potential contributions are under control in their experiments (23l . [ssj . 

B. Thermal expansion 



Even though the temperature of the cantilever clamping system is actively stabilized by a Peltier cooler, the rest 
of the system is still subject to global thermal variation as large as 1 K during the measurement runs. In order to 
understand the stability issues of our apparatus, we have studied the cantilever frequency for a time interval much 
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FIG. 14: Minimizing voltage versus sphere-cantilever distance resulting for the data analysis of various experiments. Shown 
are data from run 1 and run 2 of our experiment (top-left, with a zero minimizing potential at the closest approach for run 1, 
and a non-zero value of Vb — -150 mV for the closest approach of run 2), data from [l3| (top-right, courtesy of F. Capasso and 
H.B. Chan), results from the experiment of the Grenoble group (bottom-left, courtesy of J. Chevrier and G. Jourdan) , and 
the experiment at lUPUI (bottom-right, courtesy of R. Decca). The first three sets of data evidence a distance dependence of 
Vb upon the sphere-plane separation, which could affect significantly the accuracy without taking into account this systematic 
effect, while this behaviour is absent in the bottom-right plot. 



longer (12 hours) than the typical duration of a run (of 3-4 hours) and we have show the temporal evolution of the 
curvature coefficient Kei at a given nominal distance in Fig. 8. 

To understand the impact of thermal expansions let us consider an aluminum slab of L =1 cm with its coefficient 
of linear expansion a = 24 x 10~^/K, situated underneath the cantilever. For a typical temperature change over a 
period of 12 hours of AT — IK, the incremental expansion due to temperature is then AL = aLAT =240 nm, which 
is almost identical as the total distance variation shown in Fig. 8. This translates into a range of 5-10 nm of average 
variation in the shot- to-shot electrostatic measurement, i.e. the change in the absolute gap separation during the 
acquisition of an individual data point. To evaluate the robustness of the fits, let us assume that the last position 
registered by the PZT right before hard contact is misplaced by 8 nm either closer or away from the surface, for 
instance due to thermal expansion. One underlying assumption held throughout our analysis has been the distance 
measurement is solely controlled by the action of the piezoelectric transducer, but it is evident at this point that 
the actual distances could be also affected by thermal expansion. In order to see how this small variation due to 
temperature change influences the final results of the fitting procedure, we have adjusted the last data point of the 
closest approach by 8 nm forward or backward in an attempt to mimic the effect of a thermal expansion hypothetically 
affecting only the last collected point. 

The relocation of the single data point of the final distance has modified the entire set of fitting parameters as 
shown in Table 2. For the case of the calibration factor a, variations due to the modification of the single last data 
point were 14%, 6.4%, 18%, 21% for Run 1-4, respectively. Even larger changes in the values of the distance offset xq 
were found, up to 40 %. 

The exponent e of the electrostatic power law, however, remains stable over the forward and backward displacements 
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FIG. 15; Time behaviour of Vb in the same long run of Fig. 8. The apparent lack of correlation between this drift and the 
one presented in Fig. 8, which is mainly due to the drift in the sphere-plane gap distance, seems to imply that there are 
genuine time-dependent drifts of the minimizing potential even if the sphere-plane distance is kept constant, consistently with 
the finding in [5b |. 



of the single data point, although its overall variation is evidently greater than the fitting uncertainties at a given 
position, demonstrating how vulnerable fitting parameters are to small changes in the absolute gap distance. It should 
be emphasized that only a single data point is readjusted in this analysis. Because the thermal expansion is present 
throughout the whole measurement at all distances, what has been shown here only represents a subset of many other 
possibilities. 

Another source of drift has been identified in the minimizing potential, as evidenced in Fig. 15. This is consistent 
with recent reports of dedicated experiments studying the long term drift of the surface potential using torsional 





Run 1 


Run 2 


Run 3 


Run 4 


F 


a(Hz^V"^m°) 


3194 ± 97 


3200 ± 158 


4265 ± 155 


1749 ± 75 




xo (nm) 


35.7 ±0.1 


53.9 ±0.3 


43.5 ±0.2 


28.7 ±0.2 




e 


-1.74 ±0.01 


-1.78 ±0.02 


-1.84 ±0.01 


-1.60 ±0.01 




xVdof 


1.6 


0.5 


2.2 


7.0 




a(Hz^V-^m°) 


2805 ± 92 


3021 ± 153 


3732 ± 144 


1415 ±69 




Xo (nm) 


29.6 ± 0.1 


49.6 ±0.2 


35.7 ±0.2 


20.0 ±0.2 




e 


-1.70 ±0.01 


-1.77 ±0.02 


-1.80 ±0.01 


-1.54 ±0.02 




xVdof 


1.0 


0.8 


1.2 


7.0 


B 


a(Hz^V"^m'=) 


2788 ± 97 


3006 ± 162 


3608 ± 149 


1451 ± 78 




a;o(nm) 


28.7 ±0.2 


48.7 ±0.3 


33.9 ±0.2 


20.0 ±0.2 




e 


-1.70 ±0.01 


-1.77 ±0.02 


-1.79 ±0.01 


-1.54 ±0.02 




xVdof 


2.8 


1.6 


1.7 


11.5 



TABLE II: Stability test of the fitting parameters for the electrostatic calibrations. The middle rows present the fitting 
parameters a and a;o for the electrostatic calibrations with a power law exponent e taken to be a free parameter, such that 
Kct (X d''. F and B in the first and third rows represent hypothetical situations in which the last data point right before hard 
contact were moved by 8 nm forward (F) and backward (B) from the original position, respectively. 
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balances [56| . The presence of a minimizing potential slowly drifting in time implies severe limitations on the possibility 
to integrate for long time intervals precision force measurements. 

C. Relationship between Casimir and AFM measurements 

For the sphere-plane geometry, Casimir force experiments have used sphere radii much larger than typical in AFM 
apparatus {R ~ 100/^m versus 10-100 nm respectively) due to different motivations (mesoscopic effective surface 
versus higher image resolution, respectively). However, the physics ruling the instrument should be the same and one 
expects a smooth transition between the two regimes, and a close comparison between the accuracies reachable in 
these two regimes. Many of the systematic effects that can potentially mimic the Casimir force have been extensively 
studied in the AFM experiments [3, HO, H^, in particular for the investi gati on of adhesion and friction surface 
forces [60l . [6ll [6^ . [63j , many years before Lamoreaux's landmark experiment |l5j | . Diirig et al. have studied metallic 
adhesion and short range forces using Scanning Tunneling Microscopy (STM) [6^ . [Gq . while a number of papers 
are devoted to the subject of radiative transfer [66j . noncontact friction [btI. and dissipative interaction between a 
resonant cantilever and various sample surfaces in the dynamic AFM mode jSa. [bsI . [69l . FtoI . FtiI . [t^ , still an active topic 
of research. It would be productive to bridge the short-range AFM measurements and Casimir force measurements 
that, after the pioneering measurements of the macroscopic van der Waals force (73j . seem to have developed with 
somewhat divergent methodologies. 

The Casimir force may be considered as the van der Waals force in the long-range, retarded regime and therefore 
the study of its nature may be seen as an extension of the general studies of the previous van der Waals force 
measurements with AFM techniques. From this perspective, the measurement of forces between macroscopic surfaces 
at small distances becomes a rich field of study, including also the possible role played by the meniscus force, double- 
layer force, capillary force, hydration force, hydrophobic force, and steric, depletion, and bridge forces, all ultimately 
of electromagnetic origin, yet different in distance scaling and strength [43|. Although not all of these forces may be 
relevant in a given experimental situation, one must carefully check the presence of each of these forces and must be 
able to distinguish them, on the top of residual electrostatic forces and Casimir forces in themselves. It is plausible 
that, especially at the level of accuracy reachable in current and future high precision measurements of the Casimir 
force at small distances, some of these forces should appear already at the level of residuals from the best fits. 

A similar argument may be also applied to the electrostatic calibrations. The AFM community has already discussed 
issues such as friction, electrification, patch charges and potentials, while these have received little attention so far 
in the Casimir framework (apart from the estimates for the effect of the patch potentials discussed in [H, HI]), in 
spite of relevant analysis appearing in the recent literature [stI. [t^. Fzst . Also, Burnham et al. [bI] have evaluated 
the effect of patch charges by means of the image method, and compared with the previously observed long-range 
interaction of van der Waals force in the surface force microscopy. In [36[ , the force fluctuations between a cantilever 
tip {R = 1/xm) and a surface (both gold-coated) are interpreted in terms of inhomogeneous electric fields due to 
the presence of atomic steps, adsorbates, hillocks, pits, and other defects. The ubiquity of these electric fields may 
significantly limit precision and accuracy of Casimir force measurements. It is therefore of great relevance to aim also 
at calibrations with other physical effects, such as through radiation pressure [1^ and hydrodynamic forces in the 
case of the Casimir-Lifshitz force in the presence of liquid media (76| , and to perform detailed comparisons among the 
various calibration techniques, especially to filter out issues specific to the electrostatic calibrations. 

VI. CONCLUSIONS 

We have discussed electrostatic calibrations in a sphere-plane setup different from the previous ones for two features, 
namely a large radius of curvature of the sphere and a microresonator with large stiffness. This combination allows 
us to extend electrostatic calibrations to small gap distances in the presence of large spherical surfaces. We have 
evidenced anomalies in the expected scaling of the electrostatic force, and a distance dependent potential minimizing 
the electrostatic contribution. In the absence of a complete control of the underlying electrostatic physics, it is 
difhcult to calibrate the apparatus at the extent to meaningfully discuss the residuals signal. By using the electrostatic 
expectations, we have found evidence of the Casimir-Lifshitz force at small separation in one out of the four runs for 
which optimal working conditions were found and a complete data analysis have been performed. While we believe 
that some systematic effects are being emphasized by using spheres with large radius of curvature and small sphere- 
plane distances, we draw two conclusions from our results which may be of more general interest. Firstly, we have 
shown that the determination at all distances of the contact potential Vq is crucial, and its uncertainty can affect 
the entire data analysis procedure. This implies that electrostatic calibrations have to be performed in the entire 
range of explored distances, rather than being limited to larger separations for which the Casimir force is expected 
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to be negligible. Secondly, we have discussed the stability of the determination of the fit parameters, like distance 
offset, minimizing potential, and effective mass, coming from the electrostatic calibrations. Various groups are now 
trying to confirm or rule out the anomalies we have observed in our experimental setup ■ Apart from pointing out 
some limitations of the sphere-plane geometry, we believe that our discussion of systematic effects and data fitting 
robustness is beneficial for more rigorous data analysis in the next generation of experiments for any geometrical 
configuration, allowing to explore non-Coulombian forces, in the spirit of an experiment able to discover unknown 



physics rather than the demonstration of an a priori known effect [31 



*Present address: Department of Physics, Yale University, 217 Prospect Street, New Haven, CT 06520-8120, USA 

We thank R.L. Johnson for skillful technical support, and Q. Wei for experimental assistance. M.B.H. acknowledges 
support from the Dartmouth Graduate Fellowship program and the NSF GAAN program. W.J.K. acknowledges 
useful discussions with S.K. Lamoreaux and A. Parsegian, while D.A.R.D. acknowledges interesting discussions with 
G. Carugno, H.B. Chan, R. Decca, S. de Man, S.K. Lamoreaux, U. Mohideen, J. Munday, and G. Ruoso. R.O. 
acknowledges partial support by the NSF through the Institute for Theoretical Atomic and Molecular Physics at 
Harvard University and the Smithsonian Astrophysical Observatory, and is grateful to F. Capasso and H.B. Chan for 
kind permission to use a set of their electrostatic calibration data from [l3| for the analysis reported in Fig. 14, to J. 
Chevrier and G. Jourdan for useful discussions and for providing data reported in Fig. 14 prior to their publication, 
and to R. Decca for useful discussions and for providing data of his experiment reported in Fig. 14. 



[1] Casimir H B G 1948 Proc. K. Ned. Akad. Wet. B 51 793 

[2] Lifshitz E M 1956 Zh. Eksp. Tear. Fiz 29 94 [Sov. Phys. JETP 2 73] 

[3] Plunien G, Miiller B, and Greiner W 1986 Phys. Rep. 134 87 

[4] Milonni P W 1994 The Quantum Vacuum (San Diego: Academic Press) 

[5] Mostepanenko V M and Trunov N N 1997 The Casimir Effect and its Applications (London: Clarendon) 

[6] Bordag M 1999 The Casimir Effect 50 Years Later (Singapore: World Scientific) 

[7] Bordag M, Mohideen U, and Mostepanenko V M 2001 Phys. Rep. 353 1 

[8] Reynaud S, Lambrecht A, Genet C, and Jaekel M T 2001 C. R. Acad. Set. Paris IV-2 1287 

[9] Milton K A 2001 The Casimir Effect: Physical Manifestations of the Zero-Point Energy (Singapore: World Scientific) 
[10] Milton K A 2004 J. Phys. A: Math. Gen. 37 R209 
[11] Lamoreaux S K 2005 Rep. Prog. Phys. 68 201 
[12] Sparnaay M J 1958 Physica 24 751 

[13] van Blokland P H G M and Oveerbeek J T G 1978 J. Chem. Soc. Faraday Trans. I 74 2637 
[14] Bressi G, Carugno G, Onofrio R, and Ruoso G 2002 Phys. Rev. Lett. 88 041804 
[15] Lamoreaux S K 1997 Phys. Rev. Lett. 78 5 

[16] Mohideen U and Roy A 1998 Phys. Rev. Lett. 81 4549; Harris B W, Chen F, and Mohideen U 2000 Phys. Rev. A 62 
052109 

[17] Chan H B, Aksyuk V A, Kleiman R N, Bishop D J, and Capasso F 2001 Science 291 1941 

[18] Decca R S, Lopez D, Fischbach E, and Krause D E 2003 Phys. Rev. Lett. 91 050402; Decca R S et al. 2005 Phys. Rev. 
Lett. 94 240401 

[19] Decca R S, Lopez D, Fischbach E, Klimchitskaya G, Krause DE, Mostepanenko VM 2007 Phys. Rev. D 75 077101; 2007 
Eur. Phys. J. C 51, 963 

[20] Derjaguin B V and Abrikosova I I 1957 Sov. Phys. JETP 3 819; Derjaguin B V 1960 Set. Am. 203 47 

[21] Blocki J, Randrup J, Swiatecki W J, and Tsang F 1977 Ann. Phys. 105 427 

[22] Krause D E, Decca R S, Lopez D, and Fischbach E 2007 Phys. Rev. Lett. 98 050403 

[23] Decca R S, Lopez D, Fischbach E, Klimchitskaya G L, Krause D E, and Mostepanenko V M 2005 Annals of Physics 318 
37 

[24] van Zwol P J, Palasantzas G, and de Hosson J Th M 2007 Appl. Phys. Lett. 91 144108; 2008 Phys. Rev. B 77 035439; 

2008 Phys. Rev. B 77 075412; van Zwol P J, Palasantzas G, van de Schootbrugge M, de Hosson J Th M, and Craig V S J 

2008 Langmmr 24 7528 

[25] Jourdan G, Lambrecht A, Comin F, and Chevrier J 2007 larXiv:0712. 17671 /1 [physics.gen-ph] 11 Dec 2007 

[26] Petrov V, Hahn J, Fetter J, Petrov M, and Tschudi T 2005 Optics Lett. 30 3138; Petrov V, Petrov M, Bryksin V, Fetter 

J, and Tschudi T 2006 Optics Lett. 31 3167; Bryksin V and Petrov M 2007 JETP Letters 86 368; Bryksin V and Petrov 

M 2008 Phys. Solid State 50 229 
[27] Ludwig T 2008 J. Phys. A 41 164025 

[28] Chan H B, Bao Y, Zou J, CireUi R A, Klemens F, Mansfield W M, and Pal C S 2008 Phys. Rev. Lett. 101 030401 
[29] Fischbach E and Talmadge C L 1999 The Search for Non-Newtonian Gravity (New York: AlP/Springer-Verlag) 
[30] Gundlach J H 2005 New J. Phys. 7 205 
[31] Onofiio R 2006 New J. Phys. 8 237 



22 



Serry F M, Walliser D, and Maclay G J 1995 J. Microelectromech. Syst. 4 193 

Buks E and Roukes M L 2001 Europhys. Lett. 24 220; 2001 Phys. Rev. B 63 033402 

Chan H B, Aksyuk V A, Kleiman R N, Bishop D J, and Capasso F 2001 Phys. Rev. Lett. 87 211801 

Kim W J, Brown-Hayes M, Dalvit D A R, Brownell J H, and Onofrio R 2008 Phys. Rev. A 78 020101(R) 

Stipe B C, Mamin H J, Stowe T D, Kenny T W, and Rugar D 2001 Phys. Rev. Lett. 87 096801 

Speake C C and Trenkel C 2003 Phys. Rev. Lett. 90 160403 

Brown-Hayes M, Dalvit D A R, Mazzitelli F D, Kim W J, and Onofrio R 2005 Phys. Rev. A 72 052102 
Dalvit D A R, Lombardo F C, MazziteUi F D, and Onofrio R 2004 Europhys. Lett. 67 517 

Albrecht T R, Grutter P, Home D, and Rugar D 1991 J. Appl. Phys. 69 668; Kim B I 2004 Rev. Set. Instr. 75 5035 
Timoshenko S, Young D H and Weaver W Jr. 1974 Vibration Problems in Engineering (New York: Wyley) 
Rugar D, Marmin H J, and Guethner P 1989 Appl. Phys. Lett. 55 2588 
Lamoreaux S K 2008 E-Print arXiv:0808.0885 ^^1 [quant-ph] 6 Aug 2008 
Garcia R and Perez R 2002 Surface Sci. Rep. 47 197 

lannuzzi D, Lisanti M, and Capasso F 2004 Proc. Nat. Ac. Sci. USA 101 4019. 
Sader J E et al. 1995 Rev. Sci. lustrum. 66 3789 
Cappella B and Dietler G 1999 Surf. Sci. Rep. 34 1 

Bressi G, Carugno G, Galvani A, Onofrio R, Ruoso G, and Veronese F 2001 Class. Quantum Grav. 18 3943 
Albrecht T R ei al. 1991 J. Appl. Phys. 69 668 
Giessibl F 2003 Rev. Mod. Phys. 75 949 
Boyer L et al. 1994 J. Phys. D 27 1504 

Decca RS, Fischbach E, Klimchitskaya GL, Krause DE, Lopez D, Mohideen U, and Mostepanenko VM 2008 e-print 
larXiv:0809.3576V l [quant-ph] 

Lifshitz E M 1956 Sov. Phys. - JETP 1 73; Lifshitz E M and Pitaevskii L P 1980 Statistical Physics, Part 2 (Oxford: 
Pergamon Press) 

Palik E D (editor) 1995 Handbook of Optical Constants of Solids (New York: Academic Press) 
Chiu H-C, Chang C-C, Castillo-Garza R, Chen F, and Mohideen U 2008 J. Phys. A 41 164022 
Pollack S E, Schlamminger S, and Gundlach J H 2008 Phys. Rev. Lett. 101 071101 
Onofrio R and Carugno G 1995 Phys. Lett. A 198 365 

Chen F, Klimchitskaya, Mohideen U., and Mostepanenko V M 2004 Phys. Rev. A 69 022117 

Binnig G, Quale C F, and Gerber Ch 1986 Phys. Rev. Lett. 56 930 

Weaver J M R and Abraham D W 1991 J. Vac. Sci. Techn. B 9 1559 

Terris B D, Stern J E, Rugar D, and Mamin H J 1989 Phys. Rev. Lett. 63 2669 

Blackman G S, Mate C M, and Philpott M R 1990 Phys. Rev. Lett. 65 2270 

Burnham N A, Cohon R J, and PoUock H M 1992 Phys. Rev. Lett. 69 144 

Diirig U, Ziiger O, and Pohl D W 1990 Phys. Rev. Lett. 65 349 

Diirig U and Ziiger O 1994 Phys. Rev. B 50 5008 

Hargreaves C M 1969 Phys. Lett. A 30 491 

Volokitin A I and Persson B N 2003 Phys. Rev. B 68 155420 

Dorofeyev I, Fuchs F, Gotsmann B, and Wenning G 1999 Phys. Rev. B 60 9069 

Gotsmann B, Seidel C, Anczykowski B, and Fuchs H 1999 Phys. Rev. B 60 11051 

Loppacher C et al. 2000 Phys. Rev. B 62 13674 

Pfeiffer O, Nony L, Bennewitz R, Baratoff A, and Meyer E 2004 Nanotechnology 15 101 
Volokitin A I, Persson B N J and Ueba H 2006 Phys. Rev. B 73 165423 

Israelachvili J N and Tabor D 1972 Proc. R. Soc. London, Series A, Math, and Phys. Sci. 331 4435 
Robertson M A et al. 2006 Class. Quant. Grav. 23 2665 
Palacios-Lidon E and Colchero J 2006 Nanotechnology 17 5491 
Munday J and Capasso F 2007 Phys. Rev. A 75 060102(R) 

de Man S, Heeck K, and lannuzzi D 2008 e-print larXiv:0809.3858t fl [quant-ph]; J. Chevrier [private communication); S.K. 
Lamoreaux [private communication); G. Palasantzas [private communication) 



